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ABSTRACT 


A method of collecting sequential samples of falling hail was 
developed and used during the summers of 1973 and 1974. The col- 
lectors were mesh screens and had collecting apertures of ©.75 e 
and 0.58 a respectively. The collected hailstones were stored in 
polyethylene bottles, plastic beakers or plastic bags. The col- 
lectors were Beeaieed from vans working in conjunction with the 
Alberta Hail Project. Samples were obtained from three storms in 
1973 and six storms in 1974, 

Several methods of analyzing the collected hailstones for 
hailstone numbers and size distributions were evaluated. The most 
commonly used technique was to count and measure manually photo- 
graphs of the céllected hailstone samples. Size distributions 
produced by this method gave results similar to those obtained 
from hailpads, and different to a greater extent from those 
obtained from photographs of the hail lying on the ground. 

The resulting size distributions could be described by several 
density functions or distribution functions. An exponential 


density function 


NUDED = DE aeen wD 


provided the best deseription of the size distribution of the hail 
collected during the two years. The function varied from year to 


year. It was 
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in 1974, Exponential density functions of flux density and ‘hail- 
stone concentration were also calculated and compared with values 
found in the literature. A eeaedrd distribution was obtained for 
the collected hail. The log-normal distribution compared well 

with other log-normal distributions found in the literature when 

the hailstones considered were restricted to those with diameters 
greater than or equal to 0.5 cm. A power-law size distribution of 


the form 


was tried. The values for “ varied from 0.15 to 7.1, with a value 
of approximately 3 for the 1974 hailstones greater than or equal 
to O,o cm only. 

Three sets of hailstone samples from three differing hail- 
storms were sclected for a time-dependent analysis. Time changes 
of size spectra, hailstone flux, hailstone concentration, hailstone 
mean and median size, and variance from an exponential distribution 
were studied. The largest hailstones fell during the time of 


greatest hailstone flux. The size spectra of the samples were 
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most nearly exponential during the times of greatest hailstone 
Pbux. 

A simple hailshaft model was developed to predict upper-air 
hailstone concentrations by extrapolating backward in time the 
hailstone concentrations at the ground calculated from the 
collected samples. The model was applied to samples collected from 
a large hailstorm of 18 August 1974, The model indicated that 
the hail may have accumulated in the lower part of the updraft, and 
that the collected hailstones possibly came from two separate 
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CHAPTER 1. INTRODUCTION 


During the summers of 1973 and 1974 the author was employed 
byebr.s. P. Lozowski ofthe) University ‘of Alberta to assist’ in 
the operation of various field experiments performed in con- 
junction with the Alberta Hail Project (ALHAP). The work pre- 
sented here resulted from the sequential hail sampling experiment 
performed during the two summers. This investigation covers 


three: main areas: 


1. the size distributions of the collected hail samples 
Ze the time-varying properties of the sequential samples 
ae modelling the hailfall using the parameters of the collected 


hail samples. 

Because the internal structure of the collected sample was 
not determined, whether or not the collected particles were 
actually hailstones also remained undetermined. For this study, 
the author applies the word "hail" to precipitation consisting 
of frozen water and disregards the differences between hail, 
frozen raindrops and graupel. This approach therefore sidesteps 
the question of whether hailstones smaller than 0.5 cm in dia~ 


meter actually occur, 


Daleeoizes Distributions 
Several attempts have been made to describe the size dis- 


tribution of hail by a mathematical relation. The most prevalent 
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description is an exponential relation of the form 


=D 
Nio) SD =N. € ; J 1) (1) 


where 


NID) ED is the number of particles in a diameter interval D-3D 
to D 


Ne is the y axis intercept of the resulting curve 
iN. is the negative slope 
D is the diameter. 


Equation 1 also applies to rainfall (Marshall and Palmer, 1948). 
Exponential distributions for Alberta hail are given by 
Douglas (1960, 1963, 1964). Wis samples were collected in wire 
mesh baskets and stored in freezers after the end of the hail- 
fall. This collection method may have resulted in the melting 
of some of the smaller hailstones in his samples. Equation ] 


in Douglas' case was 


SOE S Ne. 
Noy fos 10S 30 eS 


Ludlam and Macklin (1960) give an exponential distribution 
for hailstones from a severe hailstorm in England. The sizes 
were obtained from photographs of the hail lying on the ground 
after the passage of the storm. They give a value of 
Ne = els ie for their distribution. 

Federer and Waldvogel (1975) obtained sequential hail samples 


from a Swiss hailstorm. The hailstone sizes were taken from a 


photograph of 30 seconds of hail collected on a foam rubber pad. 
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The mean distribution for the hailstones collected from this one 


storm was found to be 


-042D 


Noysp = ee 2D 


(3) 


Log-normal distributions are also applied to hail collections. 
Barge and Isaac (1973) give log-normal distributions for a 
sample of hailstones taken from 8 Alberta hailstorms in the 
summer of 1969. Since they were mainly interested in the shape 
features of the hailstones, their sample consisted mainly of 
large hailstones. In their paper, they also presented the data 
used in Douglas! exponential distributions in a log-normal form. 
Another distribution encountered in the literature is a 
power law distribution of the form 


_ 


Gti) =kD (4) 


Cay. is the concentration of the hail in the diameter inter- 
val D-dD to 0D 
© is the diameter, 
Examples can be found in Auer (1972) for hail and graupel from 
convective cloud systems over the High Plains, and Auer and 


Marwitz (1972) for hail and graupel near a thunderstorm updraft, 
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1.2 Time-varying Properties 

Although hailfall characteristics change greatly with time, 
some authors (List et. al., 1968; Charlton and List, 1972) model 
hailstorms as steady state phenomena. Douglas (1963) recognizes 
the problem in his discussion of size distributions: "The major 
source of error is likely to lie in the assumption that the size 
distribution remains invariant throughout the period of hailfall. 
Such invariance is most unlikely..." 

Williams and Douglas (1963) report the observed differences 
in three parameters; rain/hail sequence, occurrence of largest 
hail, and continuity of hail. Because their results were de- 
rived from Alberta Hail Studies (ALHAS) hailcard reports, they 
were able to examine spatial variations as well as time variations. 

Pell (1971) discusses the change with time and space of six 
hail parameters deduced from ALHAS hail card data. His results 
apply to the entire history of a hailstorm rather than bcing 
restricted to one point, so his conclusions do not relate dir- 
ectly to this study. He concludes that "...Alberta hailstorms 
seem to be most intense in their youth, becoming more and more 
dissipated as they pass middle age". 

Ulbrich (1974) derived the two parameters, J\ and No? of an 
assumed exponential distribution along with the maximum diameter 
from Topeler sc aear spectra of falling hail. He states that the 
total hailstone concentration for all sizes decreased from 
1000 ao LowLo ie and the median diameter increased from about 


0.1 cm to about 1.0 cm during a period of observation of 8.5 
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minutes. With regard tof and Nj» he says "...the temporal 
behavior of a and {\ implies a size distribution which gradually 
broadens to hailstones of larger diameters and for which the 
number of smallest diameter stones is decreasing." 

A sequential hail sampling experiment similar in concept 
to the one reported in this work was performed by Federer and 
Waldvogel (1975) in Switzerland. They were able to detect the 
passage of several hail cells during their collecting. ‘Their 
results show that of an exponential distribution showed an 
increasing trend during the course of the storm. 

Wisner et. al. (1972) present a time-dependent model of a 
hail—bearing cloud, in which they state, '"...a steady state 
assumption is appropriate until the formation of hail in the cloud." 
The model assumes exponential raindrop and hailstone size dis- 
tributions at all times and employs extensive parameterization 
of the microphysics. One of the properties was a vertical 
profile of hail concentration as a function of time. This 
output showed that the initial hail flux was small, rose to a 
maximum, then decreased to a minimum. ‘The time of the maximum 
hail flux occurrence was influenced by melting and evaporation 
terms, but it was usually in the first half of the hailstorm 
lifetime. ‘This theoretical result can be compared with the hail 


fluxes derived from this field experiment. 
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CHAPTER 2 FIELD CPERATIONS 


The sequential sampling experiment was planned to operate 
in conjunction with a mobile photography vehicle operated by 
the University of Alberta in cooperation with the Alberta Hail 
Project (ALHAP). The sampling experiment ran during the summers 
of 1973 and 1974. An account of the hail photography experiment 
has been piven by sLozowskiet sealee( LOV or 

The truck used during the summer of 1973 was a rented vehicle 
and modifications to it were not permitted. Consequently, the 
sampling equipment was of a very makeshift nature. The catcher 
consisted of a nylon mesh suspended about 70 cm above the ground 
by a light metal tubing frame (Figure 1). The nylon mesh had 
openings of about 0.5 cm. The relatively large mesh size was 
selected to allow rain to drain through quickly. No hail was 
ever observed to pass through the net. The net tapered to a 
spout in the center with an orifice diameter of 5 cm. Since any 
hail approaching that diameter would have destroyed the net, the 
spout was considered large enough to accomodate most of the hail 
expected to be encountered. The spout funnelled into the con- 
tainers used for storage, which will be described later. The 
area of the collecting aperture was 0.75 nae The slope of the 
collecting surface was fairly flat, the surface being held in 
tension by the weight of the spout only. The resulting slack- 


ness of the netting prevented any significant bouncing of the 
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FIGURE 1. The 1973 sequential sampling equip- 

ment. In the foreground is the collapsible hail 

catcher. Below the catcher are several polyethylene 

bottles containing cooled hexane used to hold the 
collected hailstones. To the left of the picture 

is the styrofoam box which stores the sample bottles 

while in transit. The large white box siin in the 
ack of the van is an extra freezer chest cooled 

by dry ice. 
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impacting hail. This catcher was collapsible and was carried 
along one wall of the van. When a stop was made for sampling or 
photography, it was set up in a convenient spot near the truck. 
Sampling could begin within 30 seconds of stopping the truck. 
No special effort was made to level the catcher, but it was 
probably within 5 degrees of horizontal in all sampling situations. 
In 1974 the University was supplied with its own truck to 
which modes otete could be made. The hail catcher was per- 
manently mounted on the roof of the van, with the spout leading 
to the interior (Figure 2). A similar nylon mesh was used but 
with slightly larger openings (0.7 cm) and with a greater per- 
centage of the netting filled in. This material was less des- 
irable than that used the previous year with respect to the 
opening size and open area, but it was the best available at the 
time. Because of the permanent mounting of the new catcher, 
the net. had to be stretched tightly to prevent sagging on the 
roof and accumulation of hail around the spout. Consequently, 
in order to prevent the hailstones from bouncing out, the net was 
given a steeper slope. The area of the catcher aperture was 
0.58 a and it stood 1 m above the roof of the van and 2.7 m 
above the ground. The net led to a 7.5 cm opening in the roof 
and thenstota flip-flop gate which couldidirect the hail, to 
either one of two containers. This gate permitted the rapid 
changing of the sample containers at accurately controlled 
times. 


Various means were used for storing the samples. The most 
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FIGURE 2, The 1974 sequential hail sampling 

equipment. The top photograph shows the hail catcher 
permanently mounted on the roof of the van. The lower 
photograph shows the interlor equipment seen through the 
door on the right side of the van. At the top of the 
photograph is the catcher drain with two outputs. Below 
the clock is a microphone to the tape recorder. On the 
left side of the doorway is an electric freezer. The 
other equipment relates to the hail photopraphy and air 


Sampling equipment. 
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frequently used method was to collect the hail in 250 ml poly- 
ethylene bottles containing hexane cooled well below the freezing 
point of water by dry ice. This method of collection should 
very rapidly freeze any liquid water in the hailstone thereby 
preserving the crystal structure and overall mass of the hail- 
stone. This "quenching" of the liquid water distinguished the 
spongy ice of the hailstone by transforming liquid portions 

into ice with many fine hairlike bubbles (Knight and Knight, 
TOGS) ee A ae problems were encountered with this method. 

Rain water is also collected along with the hail in spite of 

the holes in the mesh of the catcher, and the rain and hail tend 
‘to freeze into a single mass of ice at the bottom of the bottle. 
In a heavy shower, the bottle would fill up in about 10 seconds, 
leaving insufficient time for putting on the lid and storing 
the bottle before installing another. The hexane evaporates 
rapidly, and in spite of tight sealing of the bottles, they 
were frequently dried out when needed. This method was varied 
by using 400 ml plastic beakers with cardboard lids supplied 

by ALHAP. These were better than the bottles because of their 
increased capacity and quick lid operation, but the lids did 
not seal very well and spillage was common, For some storms 
(notably that of 18 August 1974) plastic bags were used. They 
had a large capacity, but the spongy ice was not quenched, and 
collected rainwater tended to melt the hailstones. An attempt 
was made to separate the rain from the hail in the beakers by 


inserting squares of metal screen in the beakers at the surface 
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of the hexane. The idea was that the hail would collect on the 
screen, while the rain would sink to the bottom of the beaker 
and displace the hexane up over the hail. Unfortunately, it did 
not work that way. The rainwater adhered to the surface of the 
hailstones, and the result was a frozen mass with a piece of 
screen imbedded in it. 

In both years, the samples were stored in a styrofoam box 
cooled= with dry ice and’ with a capacity of about 021 a In 
1974 the van was also equipped with an electric freezer that 
was used if the capacity of the box was exceeded. 

Although both catchers were successful in collecting seq- 
uential Mires aeres had various disadvantages. The 1973 
system was simple and relatively foolproof, but inconvenient 
and uncomfortable to operate. The 1974 system allowed better 
timing and faster operation, but the greater complexity in- 
creased the chance of mechanical breakdown. The 1975 net 
could be seen in operation, and interesting hailstones, as wel] 
as incipient trouble could be observed. The 1974 net was hid- 
den from view and could become plugged up or torn without the 
operator immediately becoming aware of it. The 1973 system re- 
quired a stationary operating location, while the 1974 system 
could be operated with the van in motion. 

The timing of the sampling was done by wristwatch in 1973, 
while in 1974, a clock was installed near the output of the 
catcher. The timepieces were synchronized every few hours 


with the clock at ALNAP by means of radio time checks. 
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The containers were all numbered and the sampling times 
and numbers were recorded on a tape recorder. The tape record- 
er broke down a few times, with a resulting loss of the data. 
The tapes were subsequently transcribed for later reference. 


Table 1 compares the two systems, 


TABLE 1. Comparison of the two sequential sampling systems. 
1973 1974 
Ty pe Portable, Permanently 
collapsible mounted 


Diameter of 
mesh holes (mm) 5.0 GeO 


Area of collecting 
aperture (m ) O25 0.58 


Height of collecting 
aperture above ground (m) Ow, te 


Maximum diameter of collectable 
hailstone (cm) ee) 6.0 


Time required to 
switch containers (s) 5 Onl 


Minimum possible sampling 
time for one container (s) 15 ] 


Number of sampling days 6 13 


Number of samples collected 58 79 
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The sampling vehicles used during the two summers operated 
together with the mobile sampling vehicles of ALHAP. These were 
directed into areas of high radar reflectivity as observed on 
the ALHAP radars. In 1974 the van was also equipped with a 
radio for monitoring the communications of the seeding aircraft 
and their controller. The procedure on a chase was to follow 
the instructions of the mobile controller at the radar console 
until the storm was in sight. Subsequently the course of action 
was determined by the observable features of the storm such as 
the precipitation shaft, cloud structure, and hail swath. 

About 6000 miles were covered in obtaining hail samples 
during 1973 and 12,000 in 1974. The useful sequential samples 


are summarized in Table 2. 
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CHAPTER 3. METHODS OF ANALYSTS 


The analysis of the hailstones was carried out in the cold 
rooms of the Departments of Mechanical and Civil Engineering 
of the University of Alberta. The eeahanasical samples that were 
used in the analysis are listed in Appendix 5. The primary 
objective of the analysis was to measure the numbers and ex- 
ternal dimensions of the hailstones in each sample, Several 
methods were employed in performing these measurements. 

The first attempt to determine the size distributions con- 
sisted of measuring each hailstone individually. However, 
holding the hailstones between the fingers melted them sig- 
nificantly and the very large numbers to be processed made this 
method quite impractical. 

The next attempt consisted of photographing each sample 
spread out on a dark background, and counting and measuring the 
images. On the suggestion of Geoff Strong (private communication) 
the facilities of the Department of Computing Science were used 
for some of the counting. One of the students of the department, 
Norman Tsang (1974), applied a computerized scanning system to 
the counting of photographed articles. The procedure involved 
line scanning an enlarged photograph of the sample with a tele- 
vision camera, followed by digital processing of the video signal. 


The processed signal was used to reconstruct the photograph on 
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a monitor. The sensitivity of the system to different light 
levels from the photograph was adjustable, so that a good rep- 
resentation of the photograph could be obtained. The recon- 
structed display consisted of a grid of 2000 by 2000 points. 
Typically, the area photographed by the film camera was about 
0.5 by 0.5 m, with the result that the resolution of the system 
was about 0.25 mm. The computer generated a number of rectangles, 
each one being just large enough to enclose a single particle 
image. The longest dimension of each of these rectangles was 
determined by the computer and a histogram of the size distri- 
bution of these longest sides was printed out. A source of 
error was introduced by the additional step of photographic en- 
largement from the negative. The blurring of the images, due 
to the film grain, added about a 20 percent uncertainty to the 
measurement of particles with a diameter of 1] mm. 

Another major problem was the overlapping of hailstones in 
the photograph. Since the computer program could not always 
resolve an overlapped image into two separate hailstones, they 
were sometimes counted as a single larger hailstone. This prob- 
lem was compounded by small ice chips and frost, which reduced 
the contrast of the photograph. In order to minimize this latter 
difficulty, most samples to be processed by the computer were 
sifted by a sieve with a 2 mm mesh size to remove the smaller 
particles. 

This scanning system was used to count most of the samples 


collected during 1973. The exceptions were those with very few 
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hailstones which were easily counted and measured by hand, or 
those where the hailstones were frozen into a block with rain- 
water and could not be separated. These samples were analyzed 

by trying to distinguish the individual hailstones on a photograph 
of the frozen sample and subsequently counting and measuring 

the images by hand, 

There were two major problems with the computerized system. 
One is the large error introduced by the inadequate resolution 
and contrast & teane photographs. The other is the arbitrariness 
of the selection of the sensitivity level for the gray scale, 
which affected the detection of the small particles and the 
dimensions of the large hailstones. This system therefore was 
not employed for the 1974 samples. 

The 1974 samples were photographed as before, but the neg- 
atives were projected onto a screen and the sizes of the hail- 
stone images were measured by hand. The scale of magnification 
was adjusted so that the smallest particles of interest (about 2 
mm diameter) could be measured with an accuracy of 20 percent. 

The hailstones were classified in size intervals of 1 mm A 
ruler was included in all of the photographs to permit an accurate 
determination of the scale. The dimension measured was the long- 
est linear dimension of the projected image. The reason for this 
procedure is as follows. A single ellipsoidal hailstone, under 
the influence of gravity, will most likely assume a position 

such that its longest dimension is parallel to the flat horizontal 


surface upon which it lies, and so parallel to the film plane. 
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The longest dimension is thus easily measured from the photograph. 
Further, it was frequently observed that if a hailstone was 
chipped, the chipping was such that the longest dimension of 

the hailstone tended to be the most preserved one, and therefore 
a better representation of the original hailstone than any other 
eneneione This anisotropic chipping of the hailstones appeared 
to result from the growth of large ice crystals with their long 
axes parallel to each other and parallel to the long axis of the 
hailstone. eines ice will tend to fracture along the boundary 

of a crystal, chipping on the Short axes of hailstones would be 
possible. Finally, the longest dimension is a representative 
dimension for deposit and transfer calculations, according to 
List and Dussault (1967). The most serious drawback of using the 
longest dimension as a representative particle diameter is that 

it gives rise to misleading volume and mass estimates, This error 
is illustrated in Table 3. The samples were sequential samples 
from the storm of 18 August 1974 which contained from 10 to 1COQ 
hailstones each. The total mass of each sample was both mea- 
sured directly and also calculated from the photographic dis- 


tributions using the formula: 
oe cg 
D t (6) 
2 ) Eoa07 P 


where 


n = number of millimeter size intervals considered 


N(D. ) = number of hailstones in size interval De to D-1 mm 
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D = image diameter in mm 


; -3 : ; 
Pa =) 50.89 g7cm! 7) taltypical idensity “for harlstones. 

TABLE 3. Comparison of sample masses obtained by weighing 

and from measured images. 

SAMPLE MEASURED MASS (¢g) CALCULATED MASS (g) RATIO 
O880 53.3 99.8 ies Ae 
0819 1egs'6) Vie 8 1.41 
0883 4.5 6.1 1,36 


The errors are caused chiefly by the non-sphericity of the hail- 
stones which would lead to an overestimation of the sample mass 
calculated from the hailstone measurements. The presence of 
frozen rainwater in the collected samples on the other hand leads 
to an overestimation of the sample mass obtained by weighing. 
Both of these problems were extant on 18 August 1974, so the 
errors in Table 5 are probably the worst to be expected. 

Another comparison of the relation of the longest image 
dimension to measurements from other sources can be obtained by 
comparing size distributions of hailstones from the storm of 18 
August 1974 which were obtained in several different ways. Data 
for this comparison came from the hail catcher samples, from 
hailpads (Strong, 1974), and from photographs of the hail on the 
ground at the end of the storm. The measurements do not apply to 
the same hailstones, but to similar samples taken under similar 
conditions. The samples were separated in space by less than 5 


meters, and were of the same 5 minutes of hailfall. Figure 3 
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shows the cumulative distributions of the samples measured by the 
three methods. Best-fitting straight lines were derived for the 
graph of the natural logarithms of the cumulative fraction of 
hailstones against their diameter. Table 4 shows three parameters 
of these straight lines for the various sampling methods. The 
hailpad measurements have been converted from dent size to hail- 
stone size using the calibration curve of Strong (1974). 

Table 4 shows that hailpads and the hail catcher give similar 
results. Using the longest dimension of a hailstone image as a 
representative dimension does not appear to result in errors 
in the size distribution exceeding those involved in the use of 
hailpads. Photographs of’ hail on the ground show more extreme 
differences from hailpads than the catcher samples do, possibly 
because the photographs were taken after the hail had ceased 
falling and when it had already melted significantly. 

We will now compare size distributions obtained by measuring 
projected photographs of hail catcher samples with those ob- 
tained by measuring the longest dimension of the hailstones 
directly with calipers. The hailstones were collected from the 
ground near Forestburg, Alberta on 22 July 1974. Table 5 and 
Figure 4 compare the two methods of measuring the sample. The 
results show excellent agreement between the two measuring 
technigues. The sample consisted of relatively large hailstones. 
Smaller hailstones would be influenced by melting in the caliper 
measurements and the discrepancy between the two measuring 


techniques would be larger. 
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FIGURE 3. Cumulative size distributions of hailstone samples ob- 
tained in three different ways. DIMENSION is the hail- 
stone diameter derived from the dent diameter for hail- 
pads, the geometric mean of major and minor axes of 
the hailstone images for the ground photograph, and the 
longest dimension for the hail catcher, 

-1 2 

Sample Y Bie ccm) Ss 

‘ fo) 

technigue ae an 

Hailpad 2.55 roles 7.02x107° 

Ground 4.93 -3.356 0.102 

photograph 

Hail tee -2.69 3.42x107° 

catcher 

TABLE 4. Comparison of three hail size distribution sampling 


techniques used on 18 August 1974, Ye is the inter- 
cept of the best-fitting straight line with ths vi 
axis.eeMins) thei) slopes ofe theystragights line. sass gS 
the variance of the difference between the sample 
cumulative distribution and the straight line. 
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TABLE 5. Comparison of sample measurements obtained by measur- 
ing photographic images of hailstones and by measur- 
ing their longest dimensions with vernier calipers. 
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FIGURE 4. Histograms of the cumilative size distributions of 
a hail sample taken from the ground near Forestburg, 
Alberta, on 22 July 1974, The caliper-measured 
histogram is displaced to the right a distance of 1.0 
cm along the x axis. 
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Knowing the accuracy of the measurements obtained using the 
projected images, the accuracy of the measurements obtained by 
the computerized line-scanning system can now be estimated. 
Figure 5 shows two histograms of hailstone size distributions for 
the same sample. One was obtained by measuring the photographic 
images manually, and the other was generated by the computerized 
line scanning system. The same photograph was used in both cases. 
The comparison is not very encouraging, because the total numbers 
and the modal sizes do not agree. The discrepancy is largely 
attributable to the overlapping of the hailstones, resulting in 
a large number of hailstones with a diameter of around 0.7 cm 
being counted as a smaller number of hailstones with diameters 
from l.O to 1.1 cm, by the computer. The larger numbers: of 
particles in the O to 0.2 cm range counted by the scanning sys- 
tem are small ice fragments and frost particles that can be 
recognized and discounted when working manually. The number of 
hailstones counted by the scanning system was 111, while there 
were 179 counted by hand. These errors were discovered quite 
bavesinethis Work, ang couldpnot be corrected due to, the poor 
focussing and lack of contrast in the photographs of the samples. 
Consequently, the size distributions for 1973 should be viewed 
with caution. 

In some cases, the samples were too conglomerated with 
frozen rainwater to permit separation of all the hailstones. 

When this happened, as many hailstones as possible were separated 


and photographed. The mass of the separated hailstones was then 
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measured as well as the total sample mass. The final number of 


hailstones in each size category was then estimated using the 


formula: 
N = Np Ms 
d Mop (7) 
where 
N is the estimated number of hailstones of diameter 
d : : 
d-1 to d in the entire sample. 

Ny is the number of separated and photographed hail- 

P stones in the diameter interval d-1l to d 
M. is the measured mass of the total sample 
es is the measured mass of the separated and 


photographed hailstones. 

This was not an entirely satisfactory procedure, because the 
frozen conglomerate consisted largely of rainwater and the smaller 
hailstones, while the larger hailstones with less mutual contact 
area would be less likely to stick together and would therefore 
predominate in the separated hailstones. Thus, the estimated 
distribution is likely to have an excess of large hailstones, but 
no way was seen to correct for this. In some cases, the procedure 
adopted was to prepare a few cross sections of the lump, and to 
measure the sizes of the individual hailstone cross-sections. 
Although the slice may not have passed through the center of all 
the hailstones in the section, enoughembryos were recognizable 

to indicate many hailstones that had been cut closc to their 


center. Using these hailstones with recognizable embryos, 
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reasonable size distributions could be obtained. Of course 
this procedure does not give the number of hailstones in the 
entire sample, but an estimate can be made based on the sample 
mass, provided there is not too much rainwater in it. The 
samples from 9 July 1974 were those most affected by large 
amounts of rainwater, 

There are two conclusions from this chapter that are impor- 
tant for thejfollowing work. First, the hail catcher and 
subsequent storage and analysis can produce results in hail 
size distributions comparable to those from hailpads. Secondly, 
the diameters of hailstones may be represented by their longest 
dimensions without unreasonable errors. For the remainder of this 
work, the term diameter will be used to mean the longest dim- 


ension of a hailstone,. 
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CHAPTER 4, SIZE DISTRIBUTIONS 


4.1 Exponential Distributions 
One way of representing size distributions of precipitation 
particles is by means of an exponential density function of the 


=i he8) 
Noo) A D = No & SD 


form: 


(8) 
For example, Marshall and Palmer (1948) found for rainfall: 
SESS 
Nico ctiby — OnUen= dD 
EO. 4) = (9) 
iy Shs sh 
where R is the rainfall rate in mm bow 
If the density function is integrated, a cumulative dis- 
tribution function is the result: 
% Nie eee 
= ln = — e€ 
(Di) tem Walkahie a A 
8) (10) 
where 
Nb) is the number of particles with diameter 2 D 
Nin = N / NK. is the total number of particles in the sample. 
.e) 


This is also an exponential distribution. The average diameter 
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of the sample is given by: 


© 


rae N(D) | 
D = D JD Nip) SD = WN (11) 


The mean volume is given by: 


= 
Ul 


ci oe 


V : /\ (12) 


This is not the volume of a hailstone of mean diameter. The 


volume median diameter Dy is givensoy (Ulbrich etal .,;or a): 


Do ~ 
| D Nw) £0 =! DN) Sd . 


The size distribution of all the hailstones collected in 
the sequential samples of the two summers is represented by the 
graph in Figure 6. This overall distribution is not representative 
either of an entire storm or of all the storms during a single 
summer. However, it is probably a fair estimate of the 
"representative"hail spectrum for Central Alberta. A straight 
line was fitted to the measured points of the distribution by 
the method of least squares giving the following parameters for 


the cumulative exponential distribution: 
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The total number of hailstones in the sequential samples is then 


NES Ny = EG 


DO 7 6 


(14) 


The variance of the collected sample distribution from the 
cumulative exponential straight line was calculated using the 


formula: 


where 


B = Afub 


and n is the number of size categories. The denominator (n-2) 
2 
is used to obtain an unbiased estimate of s in the case where 
ly ee 
two parameters, B and & = i AD are estimated. The 


9 
resulting variance was Ss = O,0199, 


Sete wot as ee : 


ws Ate in nts ek iD aagrente 1a 


Ft) 


it i's a 4 


le Up ehtnow wed, oh 4 r 
Piialtnd @ winhew <f 
Le Rey ay 
“9. setae Go aed ga 


31 


The main discrepancy between the least squares line and the 
sample points in Figure 6 occurs for diameters greater than about 
2.5 cm. There are several possibilities for this discrepancy, 
assuming that an exponential distribution represents the "true" 
distribution of the hailstones from which the sample was drawn, 

ihe Assuming that an individual sample has a double 
exponential size distribution, the oscillations may have been 
caused by the summation of several samples with different values 
of A . This is illustrated in Figure 7. The narrow distrib- 
utions may occur for individual hailfalls and will be discussed 
later. 

26 Because of the smaller numbers of the larger hailstones, 
random fluctuations in the numbers collected would produce 
relatively larger errors than for the smaller hailstones. But 
fluctuations in the numbers collected do not appear to be ran- 
dom at diameters greater than 2.5 cm. Instead, there is almost 
a continuous wave. 

Die Almost all of the larger hailstones came from two 
sampling locations of 18 August 1974. On both of those occasions, 
the catcher net was torn, and some of the hailstones, both 
large and small, could have been lost. The collecting area 
lost due to the tears was probably less than 40 percent. 

An exponential distribution was fitted to the samples for 
each year and for each sampling location with the results shown 
in Table 6 and in Figure 8. The maximum size of hailstone 
collected in 1974 was about twice that of 1973, while the total 


number of hailstones collected during 1973 was almost three times 
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Oscillation of a cumulative exponential 
distribution caused by the superposition 
of several distributions of differing 
characteristics. 
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that collected in 1974. These differences can be attributed 
mostly to chance variations in the sampling occurrences and 
locations, rather than to climatological factors, 

The results for the samples taken at single sampling 
sites in individual storms are also shown in Table 6. The values 
of A are generally greater in individual samples than they are 
for the total samples or for the yearly samples. This indicates 
that the size spectrum tends to be narrower for a single samp- 
ling location in an individual storm than it is for all the 
samples taken together. This fact is also pointed out by Auer 
and Marwitz (1972), who say: "...the size distribution of hail 
at a point from an individual storm is relatively narrow (almost 
monodisperse)....the Marshall-Palmer type of size distribution 
should not be expected throughout the storm." The cases observed 
in this work tend to corroborate the observations of Auer and 
Marwitz. The variance of the individual point samples from an 
exponential distribution is greater than that for the entire 
sample. This suggests, in agreement with Auer and Marwitz, that 
a cumulative exponential distribution is not as representative 
of a single-storm point-size distribution as it is of a sum of 
several hail samples. The relatively large variances of 0.08 
and 0.09 for 18 August 1974 are perhaps indicative of the small 
number of hailstones collected, or of sorting by size by the 
torn catcher. The value of 0.09 for 24 June 1974 does not seem 
to be an experimental error, as that was one sampling occasion 
when everything worked well. Referring to Figure 9, it can be 


seen that the large variance results principally from a decrease 
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in the number of hailstones from the expected exponential dis- 
tribution number in the size categories from 0.3 to (.6 cm. 

One feature that may be seen from Figure 8 is the flattening 
of the 1974 distribution at the smallest sizes. This is also 
evident in the curves for the individual storm samples, shown in 
Figure 9. This feature is characteristic of many of the sample 
distributions and it is poorly described by the cumulative ex- 
ponential distribution. Reasons for the absence of smaller hail- 
stones were offered in Chapter 3. From personal observation, 
there seems to be a relative scarcity of the smaller hailstones 
at the ground. A possible explanation of the absence of small 
hailstones is melting below the freezing level of the storm. 
This hypothesis is considered in detail in Appendix 2. The 
conclusion is that melting can indeed remove some of the smaller 
hailstones from an initial exponential distribution at the 
freezing level. This also results in a change in the shape of 
the observed distribution. 

Comparing the numbers of hailstones for each sample would 
be useful for estimating the relative intensities of different 
hailfalls during the sampling period. However, comparison of 
the total numbers is not very meaningful because the sampling 
times vary. This complication can be avoided by considering the 


mean flux density I of each sample: 
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where 
N(D) is the number of hailstones of diameter D-JD 
to D collected in a sampling time ft. 
A is the collector aperture area. 


Values for A and gt are given in Chapter 2 and in Appendix 3. 

The calculated values for the flux density of each sample are 
given in Appendix 3 and will be considered in detail in Chapter 5. 
Exponential distributions were fitted to the flux density values 
using the least squares method. The parameters of these expo- 
nential distributions are given in Appendix 3. These results 
indicate a range of flux densities of over two orders of magnitude 
up to almost ies oe aie 


The mean concentration C of hail (number per cubic meter) 


is given by: 


C(D) = I(bd)/V (0) és 


where 
I(D) is the flux density calculated using Equation 16 
V.(D) is the terminal velocity of hailstones of diameter D. 
The terminal velocity of a hailstone can be estimated by equating 


the drag force and the weight of the hailstone: 


it 
D (18) 
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where 
Pé is the density of the hailstone 
5 is the acceleration of gravity 


Cy is the drag coefficient of the hailetone 


Pa Ver the: densityaof caat. 
Values used in this work were P: = 0.89 g age. 9 = 981 cm 
a i= = =o =o “ 
seq =, Cy = "O23, “and WPA r=" 1Ps05 x 10 gocm™ 2. According to 


Strong (private communication) these are "representative" values 
for Central Alberta. This equation assumes spherical homogeneous 
hailstones. The concentration distributions at the ground for 
each sample are given in Appendix 3. The results for individual 
samples will be considered in Chapter 4. 

An exponential distribution of hail concentration with size 
for a severe storm in England was found by Ludlam and Macklin 
(1959), who estimated size spectra from photographs of hail lying 
on the cround following the storm. Their equation for the 


concentration distribution in the air at ground level was: 


Now) = 406°"? 


Because the hailstones may have melted on the ground before being 
photographed, tending to reduce the proportion of smaller hail- 
stones, this equation may have a flatter slope than the true value. 
Douglas (1963) also gives an exponential size distribution 

for hailstones collected in several Alberta hailstorms. His 


sampling techniques, however, tended to emphasize the larger 
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hailstones in a sample. Consequently, the slope he derived 
Ci \peto2-93-on>" : 
=e ies cm ) may be smaller than the actual value which 


occurred in his sampling occasions. 


4.2 Log-normal Distributions 


Another distribution used to describe collections of par- 


ticles is the log-normal distribution, given by: 


u(D) 
SD) ll exp (> db (19) 


where 
P(D) is: the probability ‘that a particle in the dis- 
tribution has a diameter < D 
Us is the standard normal variable 
vy 


es 
") 
a 
Lo 
3 
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n SAS 
Wye ( es (¥:-Y) is the standard deviation of the 
Cz] ' 

y'S. 


Figure 10 shows the total sample population plotted on log- 


normal graph paper. The graph shows that 16 percent of the 
collected hailstones were smaller than 0,1] cm in diameter, 
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FIGURE 10. Log-normal size distributions of hail samples 
collected during 1973 and 1974, and by Douglas 
(1960) and by Barge and Isaac (1973). Also 
shown is the distribution of the hailstones 
from the sequential samples with diameters 
greater than O.95 cm. 
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SO percent were smaller than 0.35 cm (the median size) and 90 
percent were smaller than 1.2 cm. Corresponding to the oscill- 
ations of the curve in Figure 6 there is a "kink" in the graph 
for larger hailstones, in this case at a diameter of 2.0 cm. The 
graph seems to consist of two curved segments that approach some 
maximum diameter asymptotically. It should be emphasized that 
for P 2 99.9 percent on the graph, only ten or fewer hailstones 
are involved with the result that this part yof the graph is not 
as reliable as the rest. 

Also included in Figure 10 are log-normal distributions of 
other Alberta hailstorm samples. That due to Douglas consists 
of 67 hailstone samples (Barge and Isaac, 1973). The distri- 
bution due to Barge and Isaac (1973) was derived from 1920 
hailstones collected from & Alberta hailstorms during the summer 
of 1969. If the hailstone samples of this work are reconsidered 
allowing only hail with a diameter greater than 0.5 cm (MANOBS), 
a distribution comparable to those of Douglas, and Barge and 
Tsaac results. 

Log-normal graphs of the sample size distributions of 
individual sampling occasions are shown in Appendix 4. Most of 
the graphs resemble Figure 10 in that they appear to consist of 
two curves connected at a kink. Some information about the 
hail sizes for each sampling location is given in Table 7. 


Particles smaller than 0.1 cm appeared in the samples on all but 


pies loccagione, The fraction of particlesalcss than 03:1 cm in 


diameter usually varied between O and about one fourth of the total 


collected sample. 
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The probabilities and diameters of the cumulative size 
distributions of the samples were transformed to two new var- 


tables, X and Y, using: 


= aa 


The transformed data were used to generate least squares fits 

to the graphs of Figure 10 and Appendix 4 in order to determine 
the parameters of the log-normal distribution most closely 
approximated by the data points. Table 8 gives the results of 
the calculations. The table shows that the slopes are not 
constant. According to List et. al. (1972), a constant slope 
means a constant standard deviation of the samples. This ill- 
ustrates the lack of similarity among the samples from each 
location. As will be shown in Chapter 5, the standard deviation 
also varies during a hailfall. The two samples with the closest 
fit to a log-normal distribution (minimum variance from a log- 
normal distribution) are that of 9 July and the second sampling 
location of 18 August 1974. On these two occasions, rather 


large problems arose in the collection of the samples. On 
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9 July 1974 one quantities of rainwater were collected along 
with the hail, and on 18 August 1974 the catcher net was torn. 
Collected rainwater would melt the collected hailstones slightly, 
with the smaller hailstones the most affected. The torn catcher 
net should not have had any effect on the size distribution, 


Since all sizes should have been more or less Cua Ly lost. 


4.3 Power Law Distributions 
Some authors (Tribarne and de Pena, 1962; and Sulakvelidze, 
1967) suggest that the size distribution of hail in a storm 


obeys a power law of the form: 


Etsy ere (21) 


where 
; BD x ; : 
CCR) is the concentration per m in the diameter inter- 
val D-§P to D 
D is thes diameter .of, the bails in cm 
k and GQ are constants. 


This distribution is not a cumulative distribution like the 
previous ones, but a density function. Sulakvelidze has sug- 
gested a value of 3 fora. Some values of k and A for 


Equation 21 are given in the literature. Auer (1972) gives 
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(29) 


for hail and graupel from a convective cloud system over the 


High Plains. Auer and Marwitz (1972) give 


@ 2 GIES wee (23) 


for hail and graupel in the air near a thunderstorm updraft. 

Power law distributions were fitted to the calculated con- 
centrations using a least squares method. The results are given 
in Table 9. The concentration for each size increment has been 
divided by the total concentration, so that k 1s less than l. 
This procedure allows a comparison of the calculated variance to 
that for the other distributions. 

The results indicate that the power law does not describe 
the storm samples very well, since the variance is about twice 
as large as for a log-normal distribution. Also, the slopes 
differ greatly from the Sulakvelidze value of 3. Figure 11 
shows a log-log plot for the samples of 1974. On this graph, 


those hailstones with a diameter greater than about 1 cm match 


a distribution of the type 
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FIGURE Ll. Power law size distribution of sequential 
hailstone samples from 1974. The two lines 
indicate the best-fitting distributions for 
all 1974 hailstones, and those greater than 1.0 
cm respectively. 
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quite well. Apparently, the discrepancies from this distribution 


are caused by the smaller hailstones. Knowing this, the least 


squares were recalculated using subsets of the samples with some 


of the smaller hailstones eliminated. Table 9 shows the results 


of the calculations on these truncated samples. 


These results are more encouraging. The slopes for the 


individual storms tend to be both larger and smaller than 3. 


The large 


"averaged" sample of 1974 has the slope closest to 3, 


surgesting that & = 3 is reasonable. In the cases where the 


smaller hailstones are neglected, the variance of the sample 


from a power law distribution is also reduced, as shown in 


Table 9. 
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CHAPTER 5. TIME-VARYING PROPERTIES 


The detailed study of the time-varying properties of 
hailfalls will be based mainly on three storm sample sets. 
The storms were subjectively judged to be of different sever- 
ities, based on the criteria of hail amount, hail size, duration 
of hailfall, and whether rain was falling vith the bail eethe 
presence of rain during a hailfall may indicate that some of 
the smal] hail produced by the storm melted before reaching the 
ground. The three storms selected occurred on 9 July 1974, a 
"weak'' storm; 24 June 1974, a "moderate" storm; and 18 August 
1974, a "severe" storm. The information for other storm samples 


is included in Appendices 3 and 5. 


Jelevariation of Size Distributions with Lime 

Since hailstones of differing size generally have different 
terminal velocities, the measured values of number flux and 
concentration in the air at ground level are not the same as those 
aloft. In addition, vertical shear of the horizontal wind will 
cause horizontal sorting of hailstones according to size. Fora 
stationary storm, if all the hailstones originated at one point 
and one time, the largest hailstones would reach the ground first, 
followed by the smaller ones. If the production of the hail 
continues with no change in the initial spectrum, then eventually 


the hail size distribution between the storm and the ground will 
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become homogeneous. If this storm now moves over a point on the 
ground, the initial hail collected at this point will consist par- 
tially of small hail as well as the large hail collected in the 
case of a stationary storm. Thus, the sequential samples for a 
moving storm may have an excess of small hailstones at the first of 
the sampling period compared to the sequential samples from a 
Stationary storm. This result could be modified if different 
assumptions about the hailstone distribution within the storm are made. 

Figures 12, 15 and 14 show the hailstone size spectra at 
ground level as a function of time for each of the three storms. 
Using these, an attempt was made to detect the sorting of the 
hail by size. The hailfall of 9 July 1974 (iipure 12) at the 
sampling site consisted of two events. ‘The hailfall began with 
a fairly narrow size spectrum in the first two minutes, which 
broadened in the next two minutes. The hail had almost ceased 
by the third sample, from five to six minutes. In the fourth 
sample interval, from seven to eight minutes, the hail had 
started again, with another narrow spectrum. It then diminished 
in intensity for the final time, while the spectrum remained 
quite narrow. 

The hail concentration of 24 June 1974 (Figure 13) was 
continuously decreasing in each sampling period. The beginning 
of the hailfall was not collected, which may explain the 
absence of a period of increasing hail concentration. The first 
sample had the broadest range of sizes, as well as the greatest 
concentration, The largest hailstones, those greater than about 
1.0 cm in diameter, ceased falling the soonest, disappearing 


largely after the first sample interval. The smallest hail re- 
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mained the most constant in concentration throughout the sampling 
period. Because of the lesser requirements of updrafts and 
liquid water required for its growth, small hail can be. expected 
to have a greater area of origin. 

The hailfall of 18 August 1974 is shown in Figure 14. In 
this storm, the hailfall began with hail smaller than 1.0 cm. 
There was an abrupt drop in the concentration of the small hail 
as the larger (up to 2.0 cm) hailstones began to arrive. The 


largest hail cane during the period from 9 minutes 20 seconds to 


10 minutes, which almost coincides with the greatest concentration 


of all hail. The smallest and largest hail stopped falling 
between 10 minutes and 10 minutes 20 seconds, with the inter- 
mediate sizes continuing to the end of the shower with a 
reduced concentration. 

Of the three storm samples, the one of 18 August 1974 
showed the greatest evidence of hail size sorting by the 
horizontal wind. The storm moved quite rapidly, with a speed 
of about 10 m Ca The storm of 24 June 1974 had a similar 
velocity, but the absence in the record of the first part of 
the hailfal] makes it impossible to know whether there was an 


abundance of small hail in it also. 


aie Variation of Other Parameters with Time 
The onset time of hailfall at the collecting sites was 


estimated from the tape transcripts. 
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1. On 9 July 1974 the hail began approximately 2 minutes after 
beginning sampling. During the first 2 minutes of sampling, 
rain was falling. 

2 On 24 June 1974 the hailfall began 1.5 minutes before 
sampling was started. 

ae On 18 August 1974 the hailfall began at the same time as 
the sampling. 

The samples were not of the same duration because the 

usual practice was to collect a sample sufficient to fill most of 

a sample container. The sampling intervals were therefore 

adjusted according to the rate of precipitation. 

The errors) indicated tby theterrorvbarsrinsethe following 
graphs come from three sources. 

hie There is a possible 1 percent error in the measurement of 
the effective area of the hail catcher aperture. It is 
felt that losses due to bouncing. and through the mesh holes 
probabbye*contribuvedrat most Srpercentsenror inithes@number. 

a Timing errors were introduced by approximate switching 
times and inaccurate clocks. For short samples of a 
duration of about 15 seconds, the maximum error is about 
So percent. 

ee Prrors in measuring the sizes and numbers of hailstones 
varied with the quality of the sample. These uncertainties 
ranged from greater than 100 percent for some of the 
samples of 9 July 1974 to 5 percent for some samples from 


18 August 1974. These values were derived from repeated 


analyses of the same samples. 
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The time variation of hailfall flux density (number of 
hailstones per unit area per unit time) is shown in Figure 15, 
The following features can be seen: 

Lae The point duration of hailfall varied from 6 minutes to 

17.5 minutes. <All of these values are less than the modal 

value of about 20 minutes reported by Pell (1971). Pell's 


value came from hailcard information. Although the small 


number of hailstorms considered here is probably not asignificant 


sample, the present results suggest that one should investi- 
gate the possibility that hailfall durations are over- 
estimated on hailcards. 


5 


Pa thes hailfall occurs™in bursts. Two bursts, sat about 
minutes and 10 minutes after the beginning of sampling are 
quate-prominent on, 9 July LO7A, 

3. The number /flux of the hailfall may vary by a factor of 


ten during the time of "heavy" hailfall. Heavy hailfall 


here refers to number fluxes greater than 10 hailstones 


Figure 16 shows the mass flux of ice as a function of time. 
Results for the storm of 9 July 1974 are not included, because 
the large amounts of collected rainwater obscured the mass 
measurements. The maxima and minima occur at the same times as 
those in Figure 15, but the relative maxima and minima are not as 
extreme as those for the number flux variations. Since the size 
distribution is such that the number of hailstones decreases as 
their diameter or mass increases, the greatest contributions to 


variations in number flux will come from small hailstones, which 
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will cause smaller changes in mass flux. 

Figure 17 shows the variation of the mean and median dia- 
meters of the samples with time. The variation of the two 
parameters is similar. The mean diameter increases to a max- 
imum during the hailfall, then decreases. To see whether these 
peaks correspond to those of number flux, the two variables, 
number flux and mean diameter, were correlated and the results 


are eiven in Table 10, 


DATE OF SAMPLE R 
9 July 1974 = 


0,24 
24 June 1974 0.59 
18 August 1974 0.03 


TABLE 10. Correlation coefficient R between number flux and 
mean diameter for three sets of storm samples. 

There is no evident correlation between the number flux and mean 
hailstone diameter from these three examples. A qualitative 
comparison of the two graphs (Figures 15 and 17) reveals that 
the greatest discrepancies between the two variables occurs at 
the end of the hailfall, when the number fluxes are small. The 
firstl.owmiunuces.of haulfall on 24 June 1974 were not collected: 
so the high correlation in that case may indicate that the two 
variables correlate well during the more intense middle por- 
tions of a hailfall. 

The hypothesis that the number flux is correlated with the 


maximum hailstone diameter of a sample was also tested. The 
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ro) 


ELGURE eh, Mean and median diameter vs time. The solid 
line indicates mean diameter. The dashed 
line indicates median diameter. 
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results are given in Table ll. 
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DATE OF SAMPLE R 

9 July 1974 0.88 
24 June 1974 0.70 
18 August 1974 Ona 
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TABLE 1] Correlation coefficient R between number flux and 
maximum diameter for three sets of storm samples. 


These results indicate a favorable correlation, suggesting that 
the conditions for large hail may be related to those for large 
guanlities of? hail. 

Figure 18 shows the change in the variance of the number 
flux with time. Since the variance is defined here to be a 
measure of the deviation of the sample distribution from an 
exponential distribution, a low variance indicates that the hail- 
fall is approaching an exponential distribution. The graphs 
show a variance minimum at the times of the number flux maximum. 
If reductions in number flux are due in part to dispersion of 
hailstones by the winds below a storm cell, then periods of 
greater number flux would correspond to storm samples less dis- 
turbed from the actual apeur pation within the storm. If the 
size distribution of hail within the hailstorm is exponential, 
then the variance of the distribution at the ground would be 
smaller for greater number flux, as observed. 

Figure 19 shows the variation with time of hailstone 
concentration, given by Pquation 17. The features of the graph 


are similar to those of Figure 15, but the maxima and minima 
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5 10 1S 
TIME AFTER BEGINNING SAMPLING (minutes) 
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FIGURE 
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2 3 4 
TIME AYTER BEGTNNING SAMPLING (minutes) 


18 AUGUST 1974 


5 1 Iie 20 
TIME AFTER BEGINNING SAMPLING (minutes) 


Variance of sample distributions from an 
exponential distribution. Open ends of the 
graphs indicate a sample too small to 
compute the variance (fewer than two hail- 
stones). 
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FIGURE 19. N of the best-fitting exponential] distribution 
oO . : 
for concentration measurements. The error 
bars indicate accumulated errors in collection 


and analysis. 
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are more prominent. This is a consequence of the concentration 
being a function of the terminal velocity of the hailstones. 


-3 -J : : 
cm for his exponential 


Douglas (1964) gives No aol. 2 mM 
distribution for 67 hail samples collected from 1958 to 1964. 

For small samples, the present results show that of. can vary 
Sitner Way DY up Lo a factor of ten, 

Ulbrich (1974) derived the two parameters, /\ and No? of an 
assumed exponential distribution along with the maximum diameter 
fromjthe Doppler radar spectrum of falling hail. His use of 
a vertically pointing radar implied a sampling volume extending 
oversthe entire depth of the storm, Figure 20 shows the temporal 
changes of these parameters in a hailstorm on 10 August 1966, 

Nw and ws decrease with time, while the maximum diameter 
increases, Some of the values given for and No are con- 
siderably larger than those found in the present work and in the 
literature (Douglas, 1964; and Atlas, 1963). This discrepancy 
seems to be caused by the rain also measured in the Doppler 
spectrum. Ulbrich states that the total hailstone concentration 
for all sizes decreases from 1000 ik fo FL ae and the median 
diameter increases from about 0.1 cm to 1.0 cm during the 
period of observation from 1502:39 to 1511:00. With regard to 

KL and No? he says "...the temporal behavior of No and A 
implies a size distribution which gradually broadens to hail- 
stones of larger diameter and for which the number of smaller 
diameter stones is decreasing." 


Another time variation study of the parameters NS Anil A 


was performed by Federer and Waldvogel (1975) for a hailstorm in 
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FIGURE 20 Change with time of A and Nina hailstorm 
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as observed by Doppler radar. (from Ulbrich, 
1974) 
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Switzerland on 6 July 1973. They interpreted their results to 
show the passage of 4 different hail cells, each lasting for 2 

to 4 minutes. At the beginning of each cell, mainly large hail- 
stones were found, and smaller hailstones at the end of the 
passage of each hail cell. In that storm, the number flux and 
median diameter did not increase simultaneously toward the center 
of the hailfall, contrasted with the results of this study in 
Figures 15 and 17. Their results indicate that /\ showed an 


increasing trend during the passage of the four hail cells. 


5.2 Statistical Significance 

Joss and Waldvogel (1969) estimated the sample size re- 
quired to make a good estimate of the parameters of an expo- 
nential distribution of precipitation particles. Their study 
applied to three cases of rainfall — drizzle, widespread rain, 
and thunderstorm rain. ‘they determined the product of the sampling 
area A and the sampling time t that was required to approximate 
either the precipitation rate R or the radar reflectivity 4 to 
the desired accuracy. As an example, to find an Kk or 4 value 
with a probability of 95 percent which deviates less than 10 
percent from the mean value, assuming a widespread rain 
of 1 mma dereem the product At had to be at least 1.5 noes 

Their technique was adapted to hail, using precipitation 


rates from 1.0 to 100.0 mn bre which correspond te mass fluxes 
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Oi eon COO Term min . The worst case sample in this work 
: c Be 
1s approximately a 10 second sample over an area of 0.58 m in 

es , -1 : : 
a precipitation rate of 1 mm hr. For this sampling case, the 
exponential distribution parameters estimated from the sample 
would be within 5 percent of the mean value 95 percent of the 
time. Other samples should give even better estimates. The 
conclusion is that, according to the Joss and Waldvogel criterion, 
the sampling times used in this experiment are long enough to 


give useful estimates of distribution parameters. 
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CHAPTER 6. HATLSHAFT MODELLING 


If the concentration size spectrum of hail at the ground 
is known as a function of time, the concentration which occurred 
at upper levels can then be inferred by extrapolating backwards 
in time. The upper concentrations will be modified by horizontal 
and vertical winds, changing air density, and growth or melting 
of the hailstones, as well as by the differential motions of 
hailstones of varying sizes with respect to the air. 

A simple two-dimensional model was developed to determine 
the previous concentrations of the hail samples . The model 
used the hail concentrations at the ground calculated from the 
sequential samples using Equation 17. The upper air concentrat- 
ions can be evaluated completely only if the ground concentrations 
are known for many points. en assumptions were made to 
simplify the preter of the modifying processes mentioned 
above. 

The hailstones were assumed neither to grow nor to melt. 
This assumption is not entirely unreasonable since there can be 
.no growth below the freezing level, and, for hailstones with 


a radius greater than about 0.5 cm, melting will not be signif- 


1cant. 


70 


(jee atae ae iio bl gon hh 
redala GAlg ude «tor wi Rupees of eal ode, 
a tbes tration Gav iiiide os Oy Mar eeriasinns vse . 
“06d o0i Ae Breer s es fa peda 5 isk ibiguigea ane : row 
lf audaton *6) (ian) (HR ART Yl te (we ee ndhetn tied ‘ . 
yn shite chake ste ee 

at reosgall! 9) Sayaigyah we Cahors tpn batt oplewwes 9 (gma ae 

i Ndrads. “gal en Irate ‘tea yee a *p nnien ee ie oe 

1S ahd? batehiisie taney eft: 8 atmejan iioeres, Chat re 
standenmeas or Sonne 6 eM aia wiles oolqrem 
ebaresnauee? bkyouw -8! th Rd Atel oem baeeione ai 
) whl e008 \ mit / Pease pwravo? Pe | wer tal 
LAs PAE oO Out O1K ee oi) 3s Shetipere ae 


thee a He pate uh i! tidings. © on eu ettesaltod 
of oat neds Oils. a? visniew Vian (| ine: tes @/ oes 
Shit Prnagh stad “ol wae Aes qiliwet®, Att wpked 
Vinge td) an Wile pnities ,@> 0S aielp oot aornas 


i 


A 


The hailstones were assumed to fall at their terminal veloc— 
ity relative to the air. The terminal velocity is estimated by 
Equation 15. The air density was determined from the hydrostatic 


equation and the ideal gas equation, essuming a constant lapse 


rate: 
3-1 
Pe ee) 
ates (25) 
RZ) RT, FR 
where 


AAS is the air density at z 
is the air density at the surface 


is the air temperature at the surface 


R 
F 
e is the lapse rate 
J is the altitude above the surface 
ey is the acceleration of gravity 
R is the specific gas constant for dry air. 
The following values were used in the aden Ne 200 ,—, {ie = 6.5 
K ka, gE = 930.6 cm ee R se ory 8 Bar 10°? zg oe Using these 
values, the freezing level occurs at about 3.1 km above the 
surface. 
Let the concentration of hailstones with diameter D at an 
altitude z at a horizontal distance from the sampling site x 
be C(D,x,z,t). ‘he collected samples correspond to C(D,0,0,t), 


where tis the sample time. The model then calculates C(D,x,z,t) 
5 


for other values of x, z, and t. 
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The model begins with the hailstone concentration for the 


smallest diameter interval calculated from the last sample near 


the ground — CCD, yx, 42. stp), 
where 
dD = ° 
is 0.5 cm 
pe O 
n 
Zaks _ O 
i 
t, = the midpoint of the last minute of sampling. 


At the coliection site the atmosphere is described by 


z = the surface temperature 
P. = the surface pressure 
is = the surface air density. The upper air conditions 
> 


are then calculated by mquation 25. This knowledge of upper 

air conditions allows the calculation of the terminal velocity 
of a falling hailstone at an altitude z using mtquation 18. A 
fixed time step of 1 minute is used to extrapolate the concen- 
trations backward in time and space in a Lagrangian fashion. 

This time step was chosen to correspond to the sampling interval 
used for most of the collected samples. Vectors representing the 
movenent of a hailstone with diameter D at an altitude z bere 
computed in the following ways. In the horizontal, the nailstone 
was assumed to move with the air velocity with respect to the 


ground. Hence the horizontal speed of any hailstone is given by 


Ve CZ) - x2 (26) 


where is the wind shear. ‘The form of Cis given in Figure 21. 
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In the vertical direction, the hailstones were assumed to fall 
with their terminal velocity with respect to the air. The 
assumed updraft profile is shown in Figure 21, with the maximum 
updraft velocity approximately equal to the terminal velocity lor 
the largest collected hailstone. Then the vertical velocity of 


a hailstone with respect to the ground is given by 


: 


Noa Ne (DL) -Vu (2) | 5 


where Vf) is the updraft velocity at z. Vent 2) will be neg- 
ative for hailstones falling with respect to the ground. Then 
during a small time step gt near an altitude z, a hailstone with 


diameter D will be displaced the following distances: 


eee Ve Dyas il 
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Then if the concentration at one point C(D,x,z,t) is known, the 
concentration at another point at an earlier time, t - ft, can 


be determined by 
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VW (2) 
Vi (z+5 72) 


(29) 
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This equation was evaluated at a grid of points separated by 
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thirty LOO m intervais in the horizontal, and fifty 200 m intervals 
in the vertical. 

Lhe updraft profile and the wind shear are approximations to 
the conditions which could have prevailed on 15 August 1974 at the 
Sampling site and time. ‘The wind shear was approximated from the 
storm velocity at 10,000 m as observed by the ALHAP radar. The 
updraft profile was made symmetrical, with a maximum velocity 
approximately equal to the terminal velocity of the eer eoitertee 
hailstone. Both the wind shear and the updraft velocity profile 
would have been more complicated in the real storm. The sauole 
hailstones were grouped into size categories 0.5 cm wide and all the 
hailstones within a size category were assumed to fall with the ter- 
Minal velocity of a sphere with the average diameter for that 
category. The computed concentrations for the id minutes of model 
time are shown in Figure 22. The concentrations were computed back— 
wards in time, so that the cumulative errors are greatest for the 
first minute, and decrease as time progresses The sampling site 
is in the lower right hand corner. All times are MDT. 

The model is limited in that the concentrations calculated 
represent only hail that fell at the sampling site. Other hail- 
stones may be at the same point in Space ond time, but, if they have 
different sizes, they will fall elsewhere. 

One feature of the results is the two separate areas of con- 
centration maxima at 18:15 and 18:16. The upper region would 
consist of hailstoues that are partially supported by the updraft. 
They first move into the area depicted in Figure 22 at about the 


; ; —1 
level of the updraft maximum anda fall with a velocity of Sms 
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FIGURE 22, 

Isolines of hailstone concentration produced by the hailfall 
model above and upstream of a hail sampling site. The 
samples used in the model were collected on 18 August 1974 
beginning at 18:10 MDT. The isolines are at approximately 
logarithmic intervals and enclose areas of concentration 
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with respect to the ground. Because of the large terminal velocity 
and the small concentration, this isolated area seems to consist 
of large hailstones. As seen in Figure 14, large hail fell in the 
later stages of this hailfall. The existence of local areas of 
concentration maxima may be detected in another way. “xamining 

the concentrations at the 3 km level, it may be seen that the 
concentration at the beginning time of 16:08 is 1 n. This 
concentration increases to 3 a? at 18:12, then suddenly to 30 a? 
at 116:13..0The concentration decreases to 10 ae at 16:14, then 
increases avain to 50 ikeae 18:15. This concentration is main- 
tained until about 16:13, when the concentration begins diminishing 
to O at 16:20. These regions of local concentration maxima may 

be the result of the passage of more than one hail cell. 

Another feature is the observation that hail which falls at a 
point comes from a region of greater horizontal extent in the upper 
levels. In this example, the collected hailstones extend over a dis- 
tance of about 1 km horizontally at the 4 km level at 16:16. This 
point is relevant to comvaring hail measured on the ground with 
radar measuresients aloft. 

The model is principally limited by not having a mechanism 
for the growth or melting of hail. Another improvement might be 
to use exponential distributions as the input, rather than the 


experimental results. This would give smoother concentration 


fields at the expense of greater deviations from the real sit- 


uation. 
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CHAPTER 7. SUMMARY AND RECOMMENDATIONS 


fork Summary 

The field operations of the sequential sampling experiment 
were largely successful. Sequential hail samples were collected 
in three hailstorms in 1973 and six in 1974, indicating an 
improvement in hailstorm sampling ability partially due to 
increased experience. The best results were obtained when one 
person was able to devote most of his time to the sample collect- 
ion. Collection of rainwater was the biggest difficulty with the 
equipment, but for the most severe storms, this problem was 
recuced.s ithe, large, mesh size,of) the. catcher? nettings didemot 
appear to be a major factor in the non-capture of hail. —Record- 
ing data on tape caused trouble on a few occasions. 

The analysis techniques were less successful. The manual 
method of counting and measuring photographs of hailstones 
gives useful results, but is tedious. The computerized line 
scanning system is not much faster when all the necessary photo- 
graphic steps are included, and it suffers from large calibrat- 
ion errors. Attempts at analyzing the internal structures of 


the hailstones were handicapped by the smal] hailstone sizes, 


81 


openarartigeay i@. esclee 


tim Seyee get iqee Tr) ergy 262 14 orelia rg” pio | 


- 
bet Ov11 00 Wee Helge if ed (Ai Ofheup > (oleaeseee gl Oy 


ce sndtepibal _99eF aver tre FFel a ony aga tint 


of G4 @) isso efi dia ar iGreen wtatelead wf 


> — 


vin Gotty Ten fated, os re yar deri aft 5 ym ht angee m 


7 
t 
»-jJoaiteo oj hapetow ef ai. A). eee 7) on) ey ee at side aa aye 


i 


ods die cedend thee 18 ven sah -oa¢0 7sJammlat Be nut tanh bod tag : 
ake weldaene 4143 athe fe atten Petes 677 +o) Ju 9 

goa hth aebiten *>¢+ beh ae Tsetse Gem eytel off 
-ayeset ° .Lhed I asarqnus-ade YO fl 903961 colege a) in 
—epcionace eet é WO wien) baewke Oger: og off 


leven edt .lelrse-ce a@6l S989 eo; imnieeaed © Paylene ptt 


weietmtiad te # cy: 2904 ahh B4itere~m ben gud tate san be 
: 


mihhe Bug baw spre: 9) eeuihey =) to! ot tekew lle 


i 


nde 7. wi! Lie eotp 4ehng) sion Pee at ee 


~inrdteo op won! enotted OP Bnew’ ,hetbelars oe a5. 
“) ourtie feprete). var @rixeinnga Ie eseweehh 


and by an inadequate hailstone sectioning method. The bandsaw 
used for sectioning could be used only on hailstones larger than 
about 1 cm in diameter, and the coarse teeth produced a rough 
surface, 

The size distribution studies indicated that the hail was 
essentially exponentially distributed with size, with the degree 
of fit decreasing for smaller sizes. The dividing point for the 
accurate description by an exponential distribution is at about 
a diameter of OFS cm, corresponding to the MANOBS definition of 
hail. The exponential distribution of the total hail sample 
was given by 

NCD), 20 285 259ne: ay (352) 
where Dis in cm. For samples from the individual storms, the 
coefficient of the exponent varied from 1.08 on COO. 0o ia 

Log-normal distributions also approximate the observed size 
spectra of hail, but with deviations for both large and small 
diameters of hailstones., About 65 percent of the measured 
particles were less than 0.5 cm in diameter. The log-normal 
distributions were also approximated by straight lines, and the 
slopes of these lines showed variation from 2.28 to 7.34. 

Power law distributions best fit the expected size dis- 
CrEbDuULLoOn. Ob 

COD) See aed 
when the sample considered was limited to hailstones with 
diameters greater than about 1.0 cm. When the smaller hailstones 


ss 
were included, the exponent tended to have a value around 1.5. 
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The time-dependent study did not permit generalizations to be 
made about falling hail, because of the small number of storm 
Samples studied. However, interesting deductions can be made on the 
basis of this limited information. It was found that the great— 
est values of hailstone number flux, and mean hailstone diameter 
tended to occur in the middle of a hailfall at a point. Number 
fluxes and concentrations at the ground reached maximum values of 
about 1000 hailstones m° min! and 2600 hailstones m >, These numbers 
may show a variation of a factor of ten from one sample to the next, 
so the instantaneous values may be considerzbly higher than the 
values given for 1 minute averages. The size distribution of the 
hail is most nearly exponential during the times of greatest hail- 
stone flux. ‘The changes of the size distribution of the falling hail 
with time may be partly due to the sorting of the hail with respect 
to size by the vertical shear of the horizontal wind. 

The modelling of the upper level concentration of the hail 
Samples was limited by the simplicity of the model. It snowed that 
for the particular case modeiled, the collected hailstones may have 
come from more than one distinct storm region. It also showed that 
hail at a point on the ground may have arrived from a region of 


greater horizontal extent at upper ievels. 
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7.2 Recommendations 

The only major improvement needed for the field operation 
of the experiment is some method of separating hail from the 
rain. <A lot of effort has been devoted to this problem by the 
author and by the personnel at ALHAP, but no satisfactory simple 
solution has been found, A centrifuge may work, but it seems 
unnecessarily complicated for such a simple problem. If the 
person concerned with operating the sampling system must devote 
a large part of his time to other matters, then an automatic 
sample container switcher would be a useful feature. If the sizes 
of the hail are the only interest, then the experiment could 
be considerably simplified by making the photographs of the hail 
in the field, by photographing a dark background that has been 
left out in the hailfall a known time. Schemenaur (1974) used 
a similar method for graupel. For internal analysis, the hail- 
stones will have to be collected. Federer and Waldvogel (1975) 
give a brief description of a hailstone spectrometer used to 
collect sequential samples in Switzerland. 

For the analysis of size distributions of collected hail, 
the greatest improvement would be to correctly calibrate the 
computer line scanning system. The problems of the system that 
were mentioned in Chapter 3 could be reduced by using sample 
sizes consisting of fewer hailstones. These smaller samples 
would be representative fractions of the collected samples. 
Alternatively, the collected sample could be divided into sub- 


samples, each of which would be analyzed and then added to- 
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gether. A television camera could be used to photograph the 
hail samples and the resulting Signal fed directly into the 
computer, thereby eliminating considerable noise and processing 
time. 

In order to clarify the points raised in Chapter 5 about 
the time-dependent properties of falling hail, more sequential 
storm samples should be studied. The three cases studied were 
subjectively picked to cover a wide range of storm severities, 
but perhaps they were not representative of the full diversity 
of hailstorms. The variation with time of the correlation 
coefficients among various parameters should be studied in 
greater detail. In particular, it was suggested in Chapter 5 
that hailstone number flux and mean diameter may correlate 
well only during the most intense hailfall. Although there is 
no reason for two parameters to be constantly correlated 
throughout a hailfall,there may exist correlations for part of 
the time. 

The model of falling hail needs growth and melting 
processes included in its equations. In addition, a more 
detailed wind structure and temperature profile should be used, 
and may be available from mobile radiosonde soundings. The 
model uses the hailfall at one point only, so an estimate of the 
areal extent of the fallen hail could be used to calculate more 
realistic values of actual upper level concentrations, Alter- 


natively, a set of sequential samples could be obtained from 
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Analysis of other hailstone parameters, such as axis 
ratios, surface roughnesses, and internal structure, will 
depend largely on the availability of a suitable sample con- 
sisting of large, wéll-preserved hailstones. Measurements of the 
surface roughness and axis ratios would enable better estimates 
of drag coefficients as well as an evaluation of the factors 
used in the heat transfer equations used in the study of melting 
hailstones. The analysis of the internal structure of the hail- 
stones may make it possible to determine the growth environment 
of the hailstones, but that will depend on the availability of 


suitable hailstones and analysis equipment. 
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Various problems affected the collection of the sequential 
samples, thereby affecting the analyses. The troubles were the 
following ones. 

1. 25 August 1973. Rain was collected along with hail in 
several of the samples, causing some melting and refreezing of 
the collected hailstones. This would alter the sample distrib- 
ution leaving smaller hailstones overall, and fewer small hail- 
stones than was actually the case. 

ae 27 August 1973, Very soft, wet hail was collected, 
which was quite fragile in storage. Consequently, there was a 
lot of hailstone chipping and frost in the analyzed samples. 

The effect on the distribution would be to produce an excess of 
small hailstones. 

3. 24 June 1974. The sampling location was on the edge of the 
storm and no samples were obtained from the area of the most 
intense hailfall. Consequently, the collected samples may have 
represented hailstones smaller than average for that storm, 

4. 2esuLy Los. The timing of the samples was not recorded. 
However, since the total time of sampling was recorded as 1 
minute for four samples, each was assumed to be 15 seconds long. 
D6 GFuuly 1974. Sampling took place in the final stages of 
the hailstorm. The samples collected consisted mostly of rainwater 
which melted the hailstones to some extent, and spoiled any mass 
measurements. The resulting distribution would have been altered 
in a fashion similar to that for 23 August 1973. 


6. 18 -August 1974. The catcher net was torn by the large 


SPLOTHIAD Na? CL Opa eaters jas lig tie hind estqinnm eat! seveedis 


a ita A on 
uae beds co cs. 


7 


; up 
Oi iba Me v0 3s jaateiteee niet : ea 
ie seed Gres for FT ; Ye te 
ieiai) GIR Ar Aion Sind x12 » 
at¥id: Chima "sowe) Wher , Lieeiepieeran tas tut swt teal re 
eT rn ae aaees 
bg 2 1tiw.cgn Tie! tow (alee me cee mantanh ‘ 


se. 
h BOe. at9ns 1 sAGi tara.’ ~Rpriale #¢ aliaeve shiny tet a 
oT 


othe oie band! SARE VCANSY aa Gert QA pnts 2 


400 (snheig dv ef Hhiew 10) etl ie ene Oe bare ott 


inahvaligdt as) 


byl bGt io biw nodduneh @atiges: o¢f .giPRE areas ae ¢ 
iedn WIP Be. AGH Dd VI Soha: +iwed f aaa , wwe 


a " : a ca 
eve) 9% @8irran Gaicollyd on, 4 Thayer roe of Set) tod 


Ewe 
eves - asatly “122. = i” (Hat? failure eo@ So Patines © = ae 
amass qy” hs as bo) Hi Aen Lie oe ee int See 
— 
* ee 
i ae nates. wi PRA iS pA \ietnue age sate ee "ve 


Sti? Obr63e@ Ft! i as lomigses fs **o0 , oo tye habd Te 9 
Jerengeye, wie? at i o201q Mel has fo>ot salHt. eivie 


ey 


mee. 970. Sef isea a)’ Javiow oa es do43) 1 ied o88 bad tee 


SHiogts TY") epi biaey ani Od Shido 09! GF 


Ph Ld 1 Asi ews 1) 20h we 


: oem an Ces o69 Pin seta OG Wea: 7 
7 ° _ a 7 
aes cs 


2 n iy 24 a eal 


hail. The rips were so large that particular size categories 
were probably not preferentially lost through the holes in the 
catcher net. The collection efficiency of the catcher was 
probably reduced to about one half. Sample spectra were 
probably not significantly affected, but the flux densities 


would have been reduced, 
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APPENDIX II 


THE MELTING OF FALLING HAILSTONES 
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The problem of the melting of a hailstone as it falls 
below the freezing level has been considered by several authors 
(Mason, 1956; and List and Dussault,1957). Mason (1956) 
gives the result that a particle with a radius of less than 
about 0.4 cm would melt completely on falling from a freezing 
level at 4 km above ground level in air with a lapse rate of 
6.5 C an ee 

A comprehensive treatment of the heat transfer between a 
hailstone and its environment is given by List and Dussault 
(1967). For a hailstone, there are four sources and sinks of 
heat, described by the following expressions. 


1. HNeat source from the freezing of accreted water. 


Qe = O95 vis DiRemlecth why 


ae Heat sink to warm the accreted cloud droplets. 
Ce = -).785 nV, FE Ws D Ol. Co (ToT, ) (A2) 


Die Heat sink by evaporation, condensation and sublimation, 


Qes mae Ci.8 Bialhe oe Dless i eu) UK (A3) 


4, Rate of heat transfer between a hailstone and its 


environment by conduction and convection. 
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where 
NY, nS wl 
a ie othe 


kinematic viscosity ofvair 


collection efficiency of the hailstone for 


cloud drops 


We is the 
D is thie 
oblate 


is the 


See 


Re 

Lis is the 
I 

cw 


is the 


free liquid water content of the air 

major axis diameter of the hailstone (assumed 
spherical) 

Reynolds number 

latent heat of fusion at the deposit temperature 
fraction of accreted water which freezes 


specific heat of water averaged over the tem- 


peratune gangs Teetoge 


To is the 
Ta is the 
Gris is the 
G is the 
Dwoa is the 


(Bayt is the 


Cnt Mis*tie 
b Mec the 
rw is the 
| is the 


d a 


deposit steady state temperature 

air temperature 

latent heat of evaporation or sublimation 
roughness factor 

diffusivity of water vapor in air 

saturation vapor pressure over the hailstone 
partial pressure of water vapor in ambient air 
surface ratio of spheroid: to a sphere 
heat™transier factor’ for a jspheroid 


thermal) conductivity of “air. 


When these equations are applied to melting, a few changes 


are required. The deposit temperature T, remains constant at 


OC. Because any 


d 


water drops that are shed will be warmed by the 


air and not the hailstone, Equation A2 becomes 
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Equation Al becomes 


Ga = 0.785 v We DRe = (A6) 


where We is the free liquid water produced in the wake of 
the hailstone by its\ melting. The collection efficiency in this 
case represents the divergence of the water droplets and water 
vapor in the wake of the hailstone and is approximately equal 
Tos i 

For a steady state the sume of the heat sources and sinks 


is zero, or 


Quy * Qes * Qe =() 


(A7) 


This formulation neglects the heat required to raise the tem- 
perature of the interior of the hailstone to 0 C. When Equations 
A3, A4, AS and A6 are substituted into Equation A7, the thermo- 
dynamic state of a hailstone below the freezing level can be 


described by (List and Dussault, 1967) 
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This leads to the following equation relating the diameter and 


the auxiliary function G below the freezing level: 
3 
P re SS by (A10) 


The liquid water content produced in the wake by the melting 


hailstone is given by 


! 1.68 k( Ta ala + Gane Dwa ae (Csn = en) 
Ws a 0.785 x54545y Le D4 


(All) 


The terminal speed of a falling spherical hailstone can be 


found by Equation 18. The air temperature is given by 
tes a tee jeez (Al2) 


; -1 
where ie is the lapse rate, assumed constant at 6.5 C km , 
and T is the surface air temperature. Using the above equations, 
s 
the mass melted from a hailstone in a time interval Apa is 


given by 


7D. V ‘whe (A13) 
5M = AL 4y jt o 


The change in diameter, ANB > LS eeLVeneDy, 
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and the diameter after time dt is given by 


ie ra Daag MO (Al15) 


This incremental procedure for the determination of the diameter 
of a melting hailstone as a function of time can be repeated 
until the hailstone either reaches the ground or melts com- 
pletely. 

pquattond Avil to Al6 were solved by computer with the 
following parameters. The air density was kept constant at a 
value of Dp =) si Mylan ss iii g on. The effects of relative 
humidity or already present liquid water were ignored, which 
would increase the rate of melting. The hailstones were 
assumed to be smooth spheres with a drag coefficient of 0.5 and 
a density of 0.89 ¢g cana! 

TablevAleshows the radius of the frozen part*of*a melting 


hailstone falling from the freezing level, upon reaching the 


ground. 


Tnitial radius (cm) 


O: ble OFS0e FOSS #1684 
Freezing 1.0 0 GS exoRZ6m 1036 
veye™ 2, 0 ) ) 0.22 
(kim) 
3.0 ) ) 0 0 
4.0 0 0 ) 0 


TABLE Al Radius of the frozen part of a melting hailstone at 
the ground level after falling from the freezing 


level 
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Table Al shows that for a radius less than about 0.4 cma hail- 
stone will melt completely while falling 3 km below the freezing 
level. The effects of melting become less significant as the 
diameter increases because the increased terminal velocity  re- 
sults in less time being spent in the air, and the greater mass 
requires more latent heat to melt the ice. 

The liquid water content of the air produced by the melting 
hail is an interesting by-product of the solution of the equations. 


Vv 
The values of \y are summarized in Table A2. 


Inityval 

radius (cm) 2 Oot Oundle 1s OM OS Oe Os Comer; a On Omen. 3 
-3 

Wome ari a) 2. Di OMe D 04) 48D COO eae. EEA 

TABLE A2 Liquid water concentrations produced at the ground 


in the wakes of melting hailstones,. 

Table A2 applies only in the wake of the falling hailstones, and 
so will not directly give the total water content of the air. 
However, it does show that showers of small hail may produce 
high liquid water contents, and, hence, high rainfall rates. 

The effect of melting on a cumulative exponential dist- 
ribution is shown in Figure Al. The hail size distribution at 
the freezing level is assumed to be exponential. The change in 


the smaller diameters (less than 0.5 cm) is qualitatively in 


agreement with that observed in the individual collected hail 


samples. 
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FIGURE Al. The effect of melting on a cumulative exponential 
distribution. The top curve is an initial dis- 
tribution at the freezing level, assumed exponential. 
The middle curve is the resulting distribution if 
the freezing level is at 1.0 km, and the lower 


curve is the result if the freezing level is at 2.90 


km. 
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APPENDIX III 
INDIVIDUAL SAMPLES -- 
COLLECTION LOCATIONS AND TIMES, 

SAMPLE DESCRIPTIONS AND 


DISTRIBUTICN PARAMETERS 
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APPENDIX IV 


LCG-NORMAL DISTRIBUTIONS 


FOR INDIVIDUAL STORM SAMPLES 


Ae iEae 


"7 atGrarwi 


FIGURE A2. lLog-normal size distributions of the 
at individual sample sites. The lines correspond 
to the following sampling occasions. 
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